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PREFACE 


The  following  tasks  have  been  performed  by  JAYCOR  under  contract  number 
N00014-78-C-0613: 

A.  "The  non-linear  theory  was  developed  for  the  efficiency  enhancement 
of  the  electron  cyclotron  maser,  by  tapering  the  external  magnetic  field.. 
This  work  has  been  published  as  NRL  Memorandum  Report  3983  by  P.  Sprang! e 
and  Robert  A.  Smith,  and  is  presented  in  this  report  as  an  appendix. 

B.  "The  experimental  design  aspects  have  been  considered  for  the 
effects  of  tapered  magnetic  fields  on  the  gyrotron  efficiency..  This  work 
is  presented  in  Part  I  of  this  report. 

'C.v--An  alternate  configuration  has  been  theoretically  investigated, 
for  the  utilization  of  a  Cherenkov-type  interaction  for  the  generation  of 
radiation  by  an  electron  beam  in  a  dielectrically  loaded  waveguide.  This 

t  ■ 

work  is  presented  in  Part  II  of  this  report. 


I.  DESIGN  OF  EXPERIMENT  FOR  TESTING  EFFECTS  OF  TAPERED 


MAGNETIC  FIELDS  ON  GYROTRON  EFFICIENCY 


1.  INTRODUCTION 

The  gyrotron  has  become  a  very  promising  source  for  high  power  mm  wave 

radiation.  To  date  with  gyrotron  oscillators  (gyromonotrons)  a  power  of 

212  kW  has  been  obtained  on  a  CW  basis  at  28  GHz^  ,  and  1  MW  on  a  pulsed  basis 
2 

at  100  GHz.  Extrapolations  to  CW  Megawatt  operation  at  >_  100  GHz  appears 

■5 

possible.  For  the  primary  application  presently  foreseen  for  this  type 

of  device  (electron  cyclotron  heating  in  magnetically  confined  fusion 

devices),  efficiency  will  be  of  critical  importance.  Present  devices  appear 

to  be  limited  to  efficiencies  of  less  than  50%  even  with  the  use  of  tapered 
12  4 

cavity  walls.  ’  *  Recently  methods  have  been  proposed  for  increasing  the 
efficiency  to  ~70%.  These  methods  employ  the  varying  of  the  magnetic  field 
impressed  on  the  cavity  as  a  function  of  axial  position.  One  method  utilizes 
a  step  in  the  field, 6  the  other  a  linear  taper.6  The  methods  are  predicted 
to  produce  comparable  efficiencies,  but  the  latter  appears  to  be  more  readily 
realized  in  practice. 

A  previous  JAYCOR  Report  No.  PSD200-80-Q02-FR  discusses  some  aspects  of 
the  theory  of  efficiency  enhancement  due  to  external  field  nonuniformity.  It 
explains  the  theory  behind  the  trapping  of  beam  particles  as  well  as, 
qualitatively  the  significance  of  the  various  groups  o*  particles  being 
trapped.  This  is  in  an  effort  to  explain  the  increased  transfer  of  energy 
from  the  particles  to  the  wave.  Thus  it  provides  background  for  the  orbit 
integrator  computer  program  developed  by  A.  Drobot. 


For  the  particular  case  of  a  linear  magnetic  taper  in  the  direction  of 
particle  motion,  extensive  modeling  and  computer  studies  have  been  performed.® 
Parametric  studies  were  performed  using  values  determined  by  items  already 
in  house  at  NRL,  (i.e.,  the  Varian  electron  guns,  power  supplies,  modulator, 
and  superconducting  magnets).  The  varied  parameters  consisted  of:  particle 
velocity  ratio  (a0  s  v  /v2Q);  cavity  length;  mode  numbers  n,£,;  wave  field 
strength  E0  ;  uniform  magnetic  field  strength  B  ;  degree  of  linear  magnetic 
tapering. 


The  initial  experiment  should  be  performed  using  a  cavity  with  a  length 
to  radius  ratio  of  8.  Since  this  is  not  one  of  the  longer  cavities  the  output 
power  is  expected  to  remain  above  100  kW  peak.  Power  output  of  the  device 

scales  roughly  as~l/L3  for  the  non-magnetically  tapered  device,  and  as 

2  5  3 

~ 1/L  for  the  magnetically  tapered  gyrotron. 


From  these  studies,  an  optimum  design  was  determined: 
TEq-ji  mode 


0.526  cm  cavity  wall  radius 
4.208  cm  cavity  length 
35.0  GHz  freq.  of  operation 
14.8%  AB  of  magnetic  taper 
12.6  kG  D.C.  magnetic  field  value. 


Since  the  product  of  Q  x  Pj-,eam  is  proportional  to  2ur  dr  d£  it  is  a 

measure  of  E^.  One  of  the  outputs  of  the  computer  code  is  the  value  of  this 
product  at  the  point  of  optimum  efficiency.  The  computed  value  for  the 
experiment  described  above  is  QP£ggm  =  1 26  MW  for  an  electronic  efficiency 
of  approximately  49%  with  magnetic  tapering.  The  predicted  efficiency  without 
tapering  is  only  computed  to  be  ~33.5%. 
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The  minimum,  or  diffraction  Q'7  of  this  cavity  is  approximately  given 
by  the  formula  |^(L/A0)2,  which  for  this  case  is  ~ 300.  It  is  expected 
that  an  experimentally  feasible  Q  will  be  in  the  range  of  2Qdi^. 

Two  different  methods  for  tapering  the  magnetic  field  should  be  utilized, 
depending  on  the  length  of  the  experimental  cavity  and  the  amount  of  tapering 
necessary.  The  first  method  is  to  position  a  shaped  steel  piece  around  the 
cavity  to  act  as  a  flux  gatherer.  This,  of  course,  allows  no  way  of 
actively  varying  the  amount  of  tapering  other  than  replacing  the  steel  cones. 
The  second  method  consists  of  a  coil  wound  around  the  external  surface  of 
the  gyrotron  in  the  cavity  region.  The  advantage  of  this  over  the  steel 
cone  is  that  the  current  through  the  coil  could  be  varied  to  produce  the 
exact  amount  of  tapering  to  optimize  the  efficiencies  for  certain  operating 
conditions.  For  high  degrees  of  tapering  heating  of  the  coil  windings  may 
become  a  problem;  pulsing  the  coil  may  be  required. 
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2.  DESIGN  CONSIDERATIONS 

The  mechanical  design  of  the  device  is  shown  in  Figure  1.  The  design 
permits  rapid  exchange  of  the  cavity,  output  loading,  and  input  beam  guide 
sections,  in  that  these  parts  can  be  inserted  from  the  electron  gun  end  of  the 
tube  without  the  breaking  of  welds.  We  note  that  the  components  are  positioned 
together  only  by  slight  pressure  from  the  end  flange  gasket  -  a  method 
feasible  only  because  the  electromagnetic  mode  of  operation,  TEg^ ,  requires 
no  axial  currents. 

(a)  Electron  Gun. 

8  9 

The  electron  accelerator  is  of  NRL/Varian  design  ’  ,  and  manufactured 
by  Varian.  The  design  calls  for  an  operating  voltage  of  70  kv,  with  a  beam 
current  of  8-10A.  Under  these  conditions,  the  ratio  of  perpendicular  to 
parallel  electron  velocities  should  be  1.5,  and  the  spread  in  the  parallel 
and  perpendicular  velocities  should  be  10%  and  7%, respecti vely.  Estimates 
indicate  that  these  spreads  should  not  significantly  affect  the  maser  mechanism.^0 

(b)  Beam  Input  Guide. 

For  initial  experimentation  the  input  guide  used  should  be  very 
similar  to  that  used  in  the  initial  gyrotron.  This  has  an  I.D.  near  the  gun 
of  1.067  cm  and  tapers  toward  the  cavity  to  .940  cm.  It  is  slotted  with 
azimuthal  slots  to  help  prevent  axial  wall  currents,  and  is  backed  by 
poco  graphite  as  an  absorber. 

There  is  some  question  as  to  whether  at  higher  beam  currents  there  may 
be  some  beam  interception  or  even  some  type  of  oscillation  occurring  in  the 
beam  guide.  Should  this  occur  a  second  type  of  beam  input  guide  might  be 
employed.  This  could  consist  of  alternating  discs  of  copper  rings  and  type 
2710  Ceralloy,  a  microwave  absorber. 
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In  order  to  thoroughly  study  the  impact  of  magnetic  tapering,  several 
different  length  cavities  should  be  used  during  the  course  of  the  experimenta¬ 
tion.  It  is  desirable  that  the  cavity  be  centered  on  the  magnetic  taper  such 
that  the  center  of  the  cavity  corresponds  to  a  point  in  the  tapered  magnetic 
field  equivalent  to  the  O.C.  magnetic  field  value.  To  accommodate  the 
different  length  cavities,  spacers  should  be  designed  to  maintain  the  proper 
cavity- field  alignment. 

(c)  Output  Waveguide  and  Collector. 

The  output  guide  was  chosen  to  be  a  standard  circular  waveguide 
size,  with  a  diameter  of  1.603  cm.  This  guide  has  approximately  fifty  0.159 
cm  diameter  holes  drilled  in  it  for  pumping.  It  also  serves  as  a  rudimentary 
collector  allowing  pulse  rates  of  ~10  pps. 

(d)  Output  Window. 

The  output  window  used  on  the  tapered  magnetic  field  gyrotron  should 

be  a  one  half  wave  length  resonant  window  fabricated  from  beryl lia  (Figure 

4 

2).  The  window  is  of  the  same  diameter  as  the  1.603  cm  waveguide.  Thus 
the  entire  device  from  the  output  waveguide  to  the  mode  converted  will  be 
of  uniform  diameter.  Test  results  for  th'S  window  are  shown  in  Figure  3. 

BeO  WINDOW 


Figure  2 
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( e )  Conv erter  and  Filter. 

The  gyrotron  discussed  in  this  report  was  designed  to  produce 
radiation  in  the  TEQ-|(circular)  mode.  This  mode,  while  being  very  low 
loss,  cannot  be  focused  and  is  difficult  to  handle.  It  is  therefore 
desirable  to  convert  the  radiation  into  a  rectangular  mode,  (TE^ g )  which 
can  be  focused  and  can  be  handled  using  readily  available  components.  The 
converter  chosen  was  made  by  Hitachi  Corporation  and  has  been  used  on  the 
previous  model  gyrotron.  A  mode  filter  which  absorbs  all  non  TE°  modes  will 
be  used  to  reduce  reflections. 
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II.  ANALYSIS  OF  A  BROADBAND  CHERENKOV  AMPLIFIER 


INTRODUCTION 


In  the  past  several  years,  there  has  been  widespread  interest 
in  a  variety  of  wave  generation  devices  based  on  un-neutral i zed  beams 
of  relativistic  or  mildly  relativistic  electrons.  Such  devices  are 
potentially  capable  of  producing  high  power,  with  high  efficiency,  in 
the  centimeter,  millimeter,  and  sub-millimeter  spectral  ranges.  An 
example  is  the  electron  cyclotron  maser,  which  has  produced  100  kW 
at  centimeter  wavelengths,  with  efficiencies  ranging  from  30%  -  50%. 
The  gyrotron  mechanism  involves  azimuthal  phase  bunching  of  the 
electron  beam,  which  occurs  during  resonant  interaction  of  the  Doppler 
shifted  beam  cyclotron  mode  with  a  waveguide  or  cavity  mode  (Fig.  1). 
The  bandwidth  of  the  gyrotron  mechanism  is  typically  of  order  1-2%. 

For  various  military  and  communications  applications,  it  could 
be  of  great  interest  to  develop  free-electron  amplifiers  with  wide 
operating  bandwidths  Auj  ~  <u.  In  this  report,  we  present  a  linear 
analysis  of  a  new  principle  for  wideband  amplification.  The  basic 
configuration  is  the  use  of  an  annular  electron  beam  drifting  in  a 
dielectrically-loaded  waveguide  (Fig.  2).  In  such  a  waveguide,  the 
normal  guide  modes  have  phase  velocities  asymptotic  to  c/^eT  where 
c  is  the  speed  of  light  and  e>1  is  the  dielectric  constant  (assumed 
to  be  isotropic)  of  the  waveguide  lining.  The  system  is  assumed  to  be 
immersed  in  an  axial  magnetic  field,  $  ,  sufficiently  strong  to  render 
the  motion  of  the  beam  electrons  one-dimensional.  In  this  case. 


the  amplification  mechanism  depends  on  axial  charge  bunching  in  the 
Cerenkov  interaction  of  a  longitudinal  space-charge  wave  of  the  beam 
with  the  Ez  -  componant  of  a  TM  guide  mode. 

The  interaction  we  envision  is  schematically  represented  in 
the  dispersion  diagrams  of  Fig.  3,  where 


w+ 


are  the  positive  and  negative  energy  beam  modes  and 


(2) 


are  the  waveguide  modes  for  the  case  of  a  constant  isotropic  dielectric 
material  filling  the  waveguide.  As  can  be  seen  in  Fib.  3a,  where  e=l , 
the  beam  and  guide  modes  do  not  cross  and  therefore  four  distinct  pro¬ 
pagating  modes  are  expected  on  the  entire  k  spectrum,  with  their  dispersion 
characteristics  given  by  slight  modification  of  Eqs  (1)  and  (2)  due  to 
coupling.  On  the  other  hand  however.  Fig.  3b,  when  e  f  1  and  in  particular 
when  e>  (c/V)  the  dispersion  curves  do  cross  and  therefore  growing 
modes  are  expected  due  to  the  interaction.  Most  interesting  appears  to 
be  the  case  when  e  slightly  exceeds  the  value  (c/V)  ,  since  then  strong 
interaction  is  expected  over  a  broad  band  of  frequencies. 

The  above  discussion  applies  also  to  the  topologically  equivalent 
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case  of  a  waveguide  with  partial  dielectric  filling,  for  which  the 


dispersion  curve  is  represented  by  the  broken  line  in  Fig.  3b.  Since 
this  configuration  is  more  compatible  with  e  ~  (C/V)  ,  the  analysis 
of  this  report  is  based  on  the  realistic  geometry  of  Fig.  2. 

In  Section  2,  we  present  the  general  linear  analysis  of  this 
system,  culminating  in  the  dispersion  relation  of  Eq.(ll).  In  Section  3, 
the  dispersion  relation  is  solved  analytically  in  a  simple  limiting  case. 
In  Section  IV,  we  present  the  conclusions  of  this  analysis  and  suggestions 
for  future  work  on  this  concept. 
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2.  ANALYSIS 

A.  Physical  Model 

The  configuration  and  geometrical  parameters  are  defined  in 

Figure  2.  The  waveguide  is  a  perfectly  conducting  sylinder  of  radius 

r^.  The  guide  is  lined  with  a  lossless  dielectric  (of  isotropic 

dielectric  constant  e)  in  the  space  <  r  <_  r^.  The  electron  beam 

is  annular,  filling  the  region  r-j  «  r  <  r The  regions  0  ^  r  <  r-j , 

r«  <  r  <  r~  are  in  vacuum.  An  applied  axial  magnetic  field  B  =  B  @ 
t  *3  0  0  z 

pervades  the  waveguide.  The  un-neutral ized  electron  beam  has  purely 
axial  motion;  this  corresponds  to  the  limit  Bq  °°,  so  that  transverse 
motion  is  suppressed  and  no  cyclotron  interaction  may  occur.  We  assume 
a  cold  uniform  annular  beam,  the  equilibrium  distribution  function  of 
which  is  given  by 

fQ  (r,p2)  =  nQH  (r-r-| )  H  (r2~r)  5  (pz  -  YQm V)  (3) 

2  -1 

where  yo  =  (1-V  )  and  H(x)  is  the  Heaviside  step  function,  and 

2  2  r 

where  wK  =  4-irn  e  ly  m. 

D  o  Jo 

As  discussed  in  the  introduction,  when  V>  c//F”  so  that  the  un¬ 
perturbed  beam  and  guide  modes  intersect,  we  expect  instability  (Fig.  4). 
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B.  Dispersion  Relation 

We  seek  to  derive  the  general  dispersion  relation  for  this 
instability,  in  order  to  obtain  quantitative  results  for  the  parameters 
of  interest  for  the  design  and  operation  of  a  device  to  act  as  a  broad¬ 
band  amplifier.  Such  parameters .indicated  in  Fig.  4,  are: 


•  The  optimal  frequency,  to  ,  for  which  the  amplification 
rate  is  largest. 

•  The  corresponding  spatial  growth  rate,  6k ^ . 

•  The  cutoff  frequency,  to  ,  corresponding  to  the  modified 

waveguide  curve.  0 

•  The  upper  cutoff  frequency,  ai-j ,  of  the  unstable  range 

•  The  choice  of  values  for  e,  V,  n  etc  for  which 

and/or  6kr  maximize. 

1  m  I 


To  derive  the  dispersion  relation  for  the  coupled  beam-guide 
system  of  Figure  2,  we  consider  TM  .modes  (i.e.,  no  azimuthal  de- 

U  U 

pendence) . 

The  wave  fields  are  given  by 
t  (r,z,9,t)  =  [ez(t)  e2  +  Er(r)  er  Jeikz-1wt  t 


t  (r,z,9,t)  =  B0(r)  eg  e  lkz"1uit 


where  E  (r)  is  the  solution  of  the  wave  equation 


1  d 

2  1 

k^-e(r)  h^(w,k,r) 

r  dr 

\  dr  ) 

cZ  J 

with  Er  and  BQ  given  by: 
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(7) 


Er  (r) 


fk 


e(r)  o>2/c2-k2 


d 

dr 


Ez(r) 


H 


e 


0) 

ck 


(8) 


and  Dr  =  £Er,  =  B0/y  ,  with  y  being  the  magnetic  permittivity, 
o 

The  quantity  h  (oj,k,r)  is  the  plasma  dielectric  function 


li  (u,u,r)  =  1 


2 

47 t  e 
m  e(r) 


dpz  fo(v»pz) 
Y3  (w  -  kvz)2 


where  the  equilibrium  distribution  function  ^0(r>Pz)  given  by  Eq.3 
will  be  assumed. 


Using  Eqs.  (4)  -  (9),  the  dispersion  relation  is  obtained  by 
applying  the  boundary  conditions: 

Ez  (r4)  -  0  , 

[Ez(r)]  =  [Be(r)]  =  0  {r  =  rl’r2»r3)  ’ 

where  £  t|/(r)J  indicates  the  jump  in  quantity  ty(r)  across  the  interface 
at  r.  The  seven  conditions  (10)  yield  the  dispersion  relation. 
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II 

o 
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where 


E 

2 


e  to2/c2  -  k2 


=  1 


yr 


( 


UK 


ui  -  kv 


(12b) 

(12c) 


the  last  expression  being  the  dielectric  for  a  cold  beam,  f0~  6(PZ  -  PQ) 
and  where  J^,  Ym  (I  ,Km)  are  the  ordinary  (modified)  Benel  functions  of 
order  m.  Equation  (11)  is  the  central  result  of  the  analysis.  It  must 
be  solved  numerically,  giving  the  wavenumber  (Re  k)  and  spatial  growth 
rate  (  -  Im  k)  for  real  to. 
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3-  APPROXIMATE  SOLUTION  OF  THE  DISPERSION  RELATION 

An  approximate  analytic  solution  of  Eq.  (11)  can  be  obtained  in 
the  limit  where  we  consider  a  thin  annular  beam  at  the  dielectric  inter¬ 
face;  i.e., 

r2  ~  r2  *  0;  r2i  =  r2'rl  <<;  ri  •  (13) 

The  introduction  of  this  limit  is  done  not  just  for  the  sake  of 
the  simplification  and  tractability  it  introduces  to  the  equations,  but 
also  because  it  corresponds  to  that  of  the  most  efficient  operation.  In 
fact,  since  the  waveguide  modes  outside  the  dielectric  material  are  pro¬ 
portional  to  the  modified  Benel  function,  their  amplitude  is  maximum  near 
the  dielectric  interface,  hence,  it  is  desirable  to  place  the  beam  the 
closest  possible  to  the  dielectric  to  maximize  the  coupling  coefficient. 

Under  these  assumptions,  Eq.  (11)  can  be  simplified  to  yield  the 
dispersion  relation 


1 


Il{kor2) 

w? 


Jl(kEr2)  Yo(kEr4)  -  Y1 (kEr2}  Jo(kEr4} 


.  2 


j] 


(a)  If  either  »  1  or  k^r^  =  -  r2)  «  1 ,  we  have 


J1 ( kEr2^  Yo  (kE*V  '  Y 1 ^ kEr2 ^  *VkErV 
Jo^ker2^  Yo  ^kEr4^  *  VkEr2^  Jo(kE)V 


ctn  (k£r42) 


giving 


1  Mkor2) 


k  IJk  r  _ 
0  0  0  2 


-  f-  ctn  (kEn42)  ♦  r21h!  •  0 


(b)  If  kQr2  »  1 ,  then 


vw 

yv? 


1  -  rrV  * 1  - 

c  or2 


gi  ving 


k  “  kP  Ctn  (kEr42^  +  r21h2  =  0 
o  t 


Equation  (17)  is  not  valid  near  the  cutoff  cu  .  The  result  (18)  is 
indentical  to  that  obtained  in  planar  geometry;  this  is  the  limiting 
case  of  large  aspect  ratio. 


(c)  Near  the  beam  modes,  k^r^  <_  it/2,  and  we  have 


r  -  J  t+^r42  =  -  r21hZ 

o  E 


The  cutoff  frequency  (k  =0)  is  obtained  from 


J-,(wnr2/c)  j, 
t  /  °.S-r  =  e 2  ctn 


V“or2' 


hoTz  i *1 
V  c  r2 


For  eJsr4 2  »  r2,  implying  wQr2/c  <  1 ,  we  find 


Qj  —  ■  ; 

0  2eV 


1  +  r2 


For  e55r42  <  r2  (^0^2/c  ~  1  )  ’  we  have 


oj0  =  2.405  c/r4 


Near  the  beam  mode,  we  have  from  Eq.  (19)  the  form 


(w  -  kc8) 


where  8  =  VQ/c. 


2  kj/U) 


TkJ  +  cr42 


21  y03 


Taking  k  =  oj/Sc  +  6k,  where  |6k|  <<  |k| 


we  find 


r21wb 


(5k)2Dg('^,k)  *  pfrjl 


where 


D  (u,k) 


1  1  c  +  rr 

k  2  ITT  42 

o  E 


From  Eqs.  (27)  and  (28)  we  deduce  an  important  constraint,  viz ,  that  the 
spatial  growth  rate  and  frequency  bandwidth  are  universely  related.  Large 
growth  rates  can  be  obtained  at  the  expense  of  limited  bandwidth,  and 
vice  versa.  To  see  this,  we  note  that  |(6km)3|  maximize  for 

e  »  1;  eS2  -  2  ,  (29) 

where  we  set  y~l  (i.e.,  the  beam  is  only  weakly  relativistic).  In  this 


an 


<5kl=  y  (l  -  i  /t)  i 


r2 l^b  jr  (zB2 


r42c“  2  eS! 


zJ) 


1/3 


(31 


Thus  in  the  limit  e1  >  r^  (for  e  »  1),  we  have 


U) 


U) 


on  _  /eB2 
o  Y  eB1  -  1 


=  /T 


(32 


under  the  conditions  (29)  for  which  1 6k  I  is  maximum. 

1  m1 


On  the  other  hand,  the  bandwidth  <Sw  is  given  by  the  upper  cutoff  fre¬ 
quency  gj.|  ,  where  Im<5k  =  0,  as 


.  .  _  3 

CO)  =  w-,  -  w_  =  -g— 
1  m  2y 


Cr21^b  ttS 


r42 


(cB2  -  1) 


-5/2 


1/3 


(33 


Thus  from  (28)  and  (33)  we  find  the  very  useful  relation 


5a) _  —  o  ip  8 _ 

clm6km  ~  c'  10  eB2  -  1 

which  yields  a  maximum  for  6w  at 
£  -*■  2 ,  eB2  -*•  1 ,  y  -*•  /2  ; 


(34 


(35 


the  regime  (35)  is  seen  to  be  the  opposite  limit  to  that  of  relation 


We  shall  briefly  discuss  the  consequences  of  relaxing  a  simplifying 
assumption  we  have  used  in  the  analysis,  viz_  that  the  beam  is  taken  to  be 
cold  and  uniform  in  density.  We  have  neglected  effects  of  temperature  and 
energy  shear.  Across  the  beam,  the  wave  field  is  given  approximately 
(for  r21  «  rj)  by 


dr2 


where 


-  k2  h2  (r) 


Ez(r)  =  0  , 


(36) 


"‘w  ■ 1  - 


fQ(r,p)  dp 


(37) 


Y3(oi  -  kv)3 

If  h(r)  is  a  smooth  function,  Eq.  (36)  has  the  WkB  solution 


E_(r)  *  h  2  exp  (k  /  h(r)  dr) 


./* 


(38) 


We  may  estimate  the  effects  of  energy  spread  and  self  fields  (manifested  in 
the  potential  <$>( r ) ,  by  taking  an  equilibrium  distribution  of  the  form 


(r,  p' 

1  ■  1  V 

^  -  e$(r)  exp 

(p2/2m  -  e4>(r)  -  EQ) 2 

Eth 

L  2  Eth  J 

,  (39) 


where  EQ  is  the  nominal  beam  energy  and  E^.  is  a  measure  of  thermal  energy 
spread.  Then  the  quantity  w£/y3(u)  -  kv)2  used  for  the  analysis  of  a  cold 
beam  is  to  be  replaced  by 
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OJ 


2 

po 


fl  ek6d  3k2  Eth  +  e26d2 

VmE0  ojbQ  2m(E0  +  e$0)  ^ 


(40) 


where  quantities  subscripted  "zero"  are  evaluated  at  the  center  of  the 
beam  and  where 


e60  =  mw2b0  r2 j/8  (41) 

is  the  potential  energy  drop  between  the  edge  and  the  center  of  the  beam. 
The  energy  shear  associated  with  the  self- fields  has  two  effects: 


(i)  The  effective  beam  density  (jJbo  is  reduced: 


2 

bo 


ek6d> 

/B Tw. 

o  b 


(ii)  The  effective  energy  spread  is  increased: 


(42) 


(43) 
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Summar 


In  the  analysis  of  Sections  2  and  3  we  have  demonstrated  the 
validity  of  the  dielectric-loaded  Cerenkov  wide-band  amplifier  concept. 

We  have  obtained  approximate  relations  for  the  growth  rate  and  bandwidth, 
revealing  an  important  tradeoff  between  them,  cf.  Eq.  (34).  which  must  be 
considered  in  design  applications. 

An  important  extension  of  this  work  would  be  to  consider  a  nonlinear 
calculation  in  order  to  elucidate  the  saturation  mechanism  (which  we 

*•  expect  to  be  beam  trapping)  and  ultimate  power  levels. 

k 

1 

W 

I 


26 


FIGURE  CAPTIONS 


Figure  1  -  Dispersion  curves  for  the  waveguide  mode  and  the  electron 
cyclotron  mode  in  the  electron  cyclocron  maser. 

Figure  2  -  Geometry  of  the  Broadband  Amplifier. 

Figure  3a-  Dispersion  curves  for  the  waveguide  mode  and  the  beam  modes 
when  c  =  1 . 

Figure  3b  -  Dispersion  curves  for  the  waveguide  mode  and  the  beam  modes 
when  the  waveguide  is  filled  with  a  dielectric  material  with 
e  >  (c/V)2.  The  broken  curve  refers  to  a  partially  fill eo 
waveguide. 

Figure  4  -  Wavenumber,  Re(k),  and  spatial  growth  rate,  Im(k)  for  the 

instability  (schematically).  The  solid  curves  indicate  the 
unamplified  propagating  modes. 
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THE  NON-LINEAR  THEORY  OF  EFFICIENCY  ENHANCEMENT 
IN  THE  ELECTRON  CYCLOTRON  MASER 

I.  INTRODUCTION 

Radiation  sources  based  on  the  electron  cyclotron  maser  mechanism  are  perhaps  the  most 
efficient  devices  for  generating  high  continuous  power  in  the  millimeter  and  submillimeter 
regimes.  The  first  theoretical  work1-3  on  the  cyclotron  maser  mechanism  showed  the  relativis¬ 
tic  nature  of  the  instability  while  experimental  confirmation  came  shortly  afterwards.1  5  Since 
this  early  work,  theoretical  and  experimental2'"27  research  on  the  maser  amplifier  and  oscillator 
has  indicated  that  these  devices  could  have  great  practical  values  in  areas  ranging  from  RF  heat¬ 
ing  of  fusion  plasmas  to  new  radar  systems.  For  applications  such  as  plasma  heating,  nigh 
efficiency  is  of  prime  importance 

The  operation  of  the  electron  cyclotron  maser  oscillator  can  be  viewed  qualitatively  as  th- 
exchange  of  energy  from  a  system  of  nonlinear  oscillators  (electrons)  to  the  fields  of  a  cavity. 
We  describe  briefly  the  steady  state  operation  of  a  maser  cavity  in  which  electrons  are  entering 
and  radiation  is  emitted  at  a  uniform  rate.  The  electrons  enter  the  cavity,  gyrating  about  an 
applied  longitudinal  magnetic  field.  Since  the  electron  cyclotron  frequency  is  energy  dependent, 
because  of  relativistic  effects,  these  electrons  can  be  viewed  as  nonlinear  (or  nomsochronous) 
oscillators.  The  electron  cyclotron  frequency  decreases  as  the  electron  kinetic  energy  increases. 
We  choose  the  initial  conditions  to  be  such  that  the  cyclotron  frequency  of  the  entering  elec¬ 
trons  is  slightly  lower  than  the  fixed  Doppler  shifted  frequency  of  the  cavity  field.  As  the  ini¬ 
tially  randomly  phased  electrons  drift  through  the  cavity  some  of  them  lose,  while  others  will 
Manuscript  submitted  February  26.  19"9 


i 


SPRANGLE  AND  SMITH 


gain,  kinetic  energy,  depending  on  the  initial  value  of  their  Larmor  phase  with  respect  to  the 
electric  field.  The  electrons  which  gain  energy  undergo  a  decrease  in  rotational  frequency, 
causing  them  to  move  further  from  resonance  with  the  cavity  fields.  Those  electrons  which 
lose  kinetic  energy,  on  the  other  hand,  get  closer  to  resonance  with  the  fields  and  lose  energy 
faster  and  for  a  longer  time.  Thus  it  is  possible  to  have  a  net  decrease  in  electron  kinetic 
energy  for  a  time  after  the  electrons  enter  the  cavity.  If  the  electrons  were  to  remain  within 
the  cavity  they  would  eventually  begin  to  gain  back  a  portion  of  the  net  kinetic  energy  lost. 
However,  in  designing  a  practical  cavity,  the  transit  time  of  the  electrons  is  made  equal  to  the 
time  it  takes  the  average  electron  kinetic  energy  to  reach  a  minimum.  In  the  steady  state  the  Q 
of  the  cavity  is  chosen  such  that  the  radiated  power  just  balances  the  rate  of  kinetic  energy  lost 
by  the  electrons.  For  a  comprehensive  treatment  of  the  nonlinear  and  linear  theory  of  the 
cyclotron  maser  amplifier  see  Ref.  (17).  A  number  of  excellent  review  papers8  28-30  are  avail¬ 
able  which  discuss  the  physical  mechanism  and  applications  of  the  maser  in  more  detail. 

Theoretical  and  experimental  progress  has  been  made  in  improving  the  maser  oscillator 
efficiency  by  contouring  the  walls  and  thus  the  fields  of  the  cavity.  Experimental  efficiencies  of 
—40%  at  wavelengths  of  —9  mm  and  continuous  powers  of  —20  kW  have  been  achieved/21' 

It  is  the  purpose  of  this  paper  to  analyze  the  nonlinear  particle  dynamics  in  the  cyclotron 
maser  and  address  the  problem  of  efficiency  enhancement.  We  show  that  by  appropriately  con¬ 
touring  the  applied  magnetic  field  or  the  cavity  electric  field  as  a  function  of  axial  position,  the 
beam  electrons  can  be  prebunched  in  phase,  resulting  in  a  significant  increase  in  efficiency  over 
the  already  high  intrinsic  value.  Although  we  analyze  a  maser  oscillator  configuration,  our 
method  of  enhancing  the  efficiency  is  applicable  to  a  maser  amplifier  with  minor  changes. 
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II.  ANALYSIS 

a.  Maser  Cavity  Mode! 

We  shall  consider  a  steady  state  maser  cavity  oscillator  with  the  geometry  shown  in  Fig.  1. 
The  electron  beam  drifts  in  the  2-direction  and  gyrates  about  the  applied  axial  magnetic  field 
B0e..  The  beam  consists  of  a  thin  current  sheet,  uniform  in  the  /-direction,  with  the  guiding 
centers  of  the  electrons  initially  uniformly  distributed  on  the  x  -  0  plane.  For  a  sufficiently 
low-current  beam,  we  may  take  the  cavity  fields  to  be  given  by  the  unperturbed  normal  modes 
of  the  cavity.  For  definiteness  we  assume  that  the  beam  drives  a  single  TE0  mn  mode.  The 
only  non-vanishing  field  components  for  this  mode  are 

Ey(x,  z,  r)  -  Y.Em.n  cos(av„?  +  d>„ „ )sin ( ) sin (km (a  -  x)) 

m.n 

„  cKn 

Bx (x,  t)  -  E„„  sin  (av„r  +  (f.m„)cos(  A„:)sinU„  (a  -  x))f 

m.n  ^m.n 

ck 

B.(x,  2.  /)  «*  £  — —  Em„  sin(u>m„r  +  <P„,„)sin( Kn:)cos(km(a  -x)),  (la,b,c) 

m.n  w m.n 

where  and  <b„„  is  the  mode  amplitude  and  phase,  km  =  mir/la,  Kn  -  nir/L,  cu„„  - 
c  yjkt  +  K}  and  m,n  —  1.  2,  3 .  For  strong  coupling,  the  electric  field  must  be  max¬ 

imum  at  the  beam  location,  i.e.  m  odd,  and  have  a  frequency  close  to  the  electron  cyclotron 
frequency. 

b.  Particle  Dynamics 

The  particles  enter  the  cavity  at  r  -  0,  randomly  phased,  with  the  same  axial  velocity 
v„.  -  f3„.c  and  the  same  transverse  velocity  v,_  -  /3,  c.  In  order  to  obtain  tractable  nonlinear 
orbit  equations  we  make  a  number  of  simplifying  assumptions,  all  of  which  are  physically  well 
founded  We  assume  that  the  transverse  wavelength  of  the  cavity  mode.  2ir, /c„,  is  large  com¬ 
pared  to  the  electron  Larmor  radius  as  well  as  the  excursions  of  the  particle  guiding  centers  in 
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the  x-direction.  hence  \kmx\  «  1,  sin(/cm(a  -  x))  =  sin(/cma)  and  cos (k,„(a  -  x))  = 
xkms\n(kma)  for  m  odd.  The  particle  velocities  are  considered  to  be  mildly  relativistic,  i.e.. 
(:/c)2,  (x/c)2  and  (j)/c)2  «  1.  Finally  we  assume  that  the  cavity  is  long,  i.e.,  L  »  2a.  and 
that  the  cavity  field  amplitude  is  small  compared  to  the  applied  magnetic  field,  i.e.. 
\EmJB„\  «  1. 

For  a  single  TE  cavity  mode  the  orbit  equations,  in  the  absence  of  self-fields,  take  the 
approximate  form 

7,  {yjc)  = 

(yy)  ~  fi„x  —  £  sin(Aa)sin(3fe)cos(cor  +  d>) 

dt  m„ 

^r~  ~  0  (2a,b,c) 

at 

where  y  -  (1  -  Or2  +  y2  +  i2)/c2)"l/2  >  1,  (l0  -  \e\B0/m0c.  The  subscripts  on  Em„ ,  km,  K ,, 

o)„n  and  <t>m„  have  been  dropped.  In  obtaining  Eqs.  (2),  some  higher  order  terms  have 

been  neglected,  these  terms  are  typically  smaller  than  the  dominant  terms  by  the  factor 
(i/c)(cK /oi)(xk)  «  1. 

From  Eqs.  (2a)  and  (2c)  we  find  that  the  canonical  momenta  in  the  x  direction  as  well  as 
the  axial  velocity  of  the  particles  in  the  cavity  are  approximately  conserved,  i.e.,  yx  +  fl„.v  -  0 
and  i  —  v„..  Making  use  of  these  two  approximate  constants,  Eq.  (2b)  can  be  put  into  the 
form 

y  +  0|(1  -  /32.  -  niZ/c2  -  rVc2)^ 

—  -  £  sin(/ca)sin(Afv<,r(r  -  t0))cos(wi  +  <i>) .  (3) 

where  r„  is  the  entry  time  of  the  particle  into  the  cavity  and  we  have  used  the  approximation 

y  ~  1  +  (/32;  +  n2y2/c2  +  >',2/c2)/ 2.  The  phase  bunching  process  in  the  cyclotron  maser 
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occurs  on  a  time  scale  which  is  slow  compared  to  the  cavity  wave  period  or  electron  cyclotron 
period.  Since  Kvo;  «  cu  =  0„  we  see  that  the  solution  to  Eq  (3)  can  be  written  in  the  form 

y  —  r(i)sin(<t>t  4-  9(t))  (4) 

where  rand  9  vary  slowly  in  time,  i.e.,  (r/r|,  9  «  cu.  Substituting  Eq.  (4)  into  (3)  and  using 

the  approximations  fljy2  +  >-2  =  ca2>-2  and 

fl2(l  —  /32.  —  (i? P-/c)  —  o?  ~  —  2 cu{cu  —  ft,  [1  —  (j32..  +  (u2r2/c2)/2]), 
we  find  that  r  and  9  satisfy  the  following  equations: 

r  —  -  ec  sin (Kv0.t)cos(9  -  6), 

9  -  ft„/y(r)  -  cu  +  ~  sin(Af v„.r)sin(d  -  <t>),  (Sa,b) 

where  y(r)  -  1  +  l/2(/3|.  +  <j2r/r)  >  l  and  e  3  (£/2£0)sin(/cn)  «  1.  In  deriving  Eqs. 
(5)  we  have  neglected  small  terms  in  r,  9,  r9  and  I92.  Note  that  the  particle  time  of  entry.  t0, 
can  be  arbitrarily  set  equal  to  zero  without  loss  of  generality 

Since  fi ,  =  cu,  the  term  ojr/ccan  be  identified  as  the  magnitude  of  the  transverse  particle 
velocity.  The  kinetic  energy  of  a  particle  in  the  cavity  is  £  »  (y  -  Dm, c2  and  is  a  function  of 
r  only.  The  evolution  of  r,  given  by  the  simultaneous  solutions  of  Eqs.  (5),  gives  the  particle 
energy  as  a  function  of  time  (axial  position)  in  the  cavity.  The  nonisochronous  nature  of  the 
maser  process  is  evidenced  by  the  fact  that  the  electron  cyclotron  frequency  (the  first  term  on 
the  RHS  of  Eq.  (5b))  is  energy  dependent  because  of  the  relativistic  term  a >V/'c2.  The  energy 
dependent  cyclotron  frequency  is  responsible  for  the  phase  bunching  of  the  particle  ensemble, 
which  is  initially  distributed  randomly  in  9. 

c.  Forward  or  Backward  Wave  Coupling 

The  standing  wave  associated  with  a  cavity  mode  is  the  superposition  of  a  forward  and  a 
backward  propagating  wave.  We  may  further  simplify  the  description  of  the  nonlinear  particle 
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dynamics  contained  in  Eqs.  (5)  by  noting  that  under  certain  conditions  the  beam  particles  will 
predominantly  couple  to  either  the  forward  or  backward  wave.  This  will  permit  us  to  derive  a 
useful  constant  of  the  motion  for  analyzing  the  nonlinear  dynamics  of  the  particles.  To  see  this 
we  define  the  variable  |±  *■  0  ±  K\0.i  and  note  that  for  small  cavity  field  levels,  i.e.,  «  — •  0, 
is  appropiately  given  by 

=  Acu.  1  +  0„ 

where  Su±  =w -Cla/y(r0)  =F  K\a;.  Interaction  with  the  forward  wave  implies  |Aiu_|  » 
|Atu,|  =  0  while  for  the  backward  wave  resonance  we  would  have  |AwJ  »  |Aw_|  =  0.  If, 
however,  |A<j+|  =  |Acu_i,  both  waves  couple  to  the  beam  about  equally  and  the  overall  interac¬ 
tion  is  weak,  resulting  in  low  efficiency.  It  is,  therefore,  desirable  and  possible,  by  properly 
choosing  the  system  parameters,  to  have  one  of  the  wave  interactions  dominate  over  the  other. 
If  this  is  done  the  nonlinear  coupled  equations  in  (5a, b)  can  be  put  into  the  form 

r  -  T  -y  cos 

k±  -  O 0/y(r)-w  ±Kvo:  ±  —  sin^  (6a, b) 

where  without  loss  of  generality,  <b  has  been  set  equal  to  3ir/2  in  Eqs.  (6)  and  the  upper 
(lower)  sign  corresponds  to  forward  (backward)  wave  resonance.  These  equations  are  identical 
in  form  to  those  describing  the  particle  orbits  in  a  cavity  field  that  is  uniform  in  For 
definiteness  we  will  consider  only  the  forward  wave  resonance  case  (upper  sign);  this  implies 

that  oj  >  n„/y.  With  obvious  changes  in  parameter  definitions  our  results  apply  directly  to  the 
backward  wave  resonance  case.  Dropping  the  subscript  on  £*  we  get 


f  -  no/y(r)  -d>  +  y  sin  (7a, b) 
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where  J>  -  iu-Kvo:.  This  system  of  equations  is  equivalent  to  that  obtained  from  the  well 
known  Duffing  equation  on  a  slow  time  scale."  The  nonlinear  orbits  defined  by  Eqs  (7)  can  be 
shown  to  have  the  following  constant  of  motion 

C(r.f)  -  (d,-ajy<r))2  -  «  —  sin  {  (8) 

n„u>  c 

Typical  curves  in  (r,f)  space  are  shown  in  Fig.  (2)  for  various  values  of  C.  In  general,  particles 
may  lie  on  open  or  closed  (trapped)  orbits. 

Particles  starting  on  a  particular  constant-of-motion  curve  remain  on  this  curve.  The  rate 
at  which  the  particles  move  along  the  C-curves  is  governed  by  the  solutions  of  the  orbit  equa¬ 
tions  in  (7).  We  assume  no  interaction  among  the  particles  themselves  (/.«?.,  we  neglect  spac^- 
charge  effects).  Thus  we  may  represent  the  beam  by  an  ensemble  of  a  small  number  of  parti¬ 
cles.  typically  32.  Nunu'ically  solving  Eqs.  (7).  for  such  an  initial  ensemble  of  particles,  the 
distribution  of  particles  for  various  times  (or  axial  positions  within  the  cavity,  Is  shown  in  Fig. 
(3).  The  particles  enter  the  cavity  uniformly  distributed  in  £  over  an  interval  of  2tt  radians,  all 
with  the  same  initial  value  of  transverse  velocity,  i.e..  ,  - r, ,  see  Fig.  (3a)  As  time  progresses, 
the  particles  move  along  their  respective  C-curves.  In  the  linear  regime  (Fig.  3b)  a  significant 
amount  of  spreading  in  r  along  with  a  small  degree  of  bunching  in  £  has  developed  As  the 
particles  begin  to  bunch  in  phase  angle  £  (Fig.  3c)  most  of  them  are  contained  within  the 
energy  loss  regime  — w/2  <  f  <  tt/2  and  considerable  spreading  in  r  has  occurred.  The  parti¬ 
cles  are  distributed  over  a  wide  range  of  C-curves.  with  two  ensuing  consequences  First,  as 
the  particles  lose  energy  they  move  out  of  the  energy  loss  region,  -n  2  <  t  <  rr  2.  Secondly 
not  all  the  particles  reach  their  minimum  energy  at  the  same  time  Figure  f  3d )  shows  the  parti¬ 
cle  ensemble  in  the  saturated  state  where  the  particles  are  deeply  trapped  and  the  average 
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transverse  kinetic  energy  is  ai  a  minimum.  If  permitted,  the  particle  ensemble  would  continue 
to  rotate  and  regain  the  lost  kinetic  energy. 

The  efficiency  of  energy  extraction  from  the  particles  is  defined  as 

V<l)  -  N;'  £  (y(r,  (f))-y(r„))/(y(r0)-l)  (9) 

i-i 

where 

y(r,(i))  -  y(r„)  +  w2c~2(r2(r)  -  r,})/2 

(c.f.  below  Eqs.  (5)  ),  r0  —  r,(r— 0)  and  r,(t)  is  the  Larmor  radius  of  the  /  th  particle  which 
entered  the  cavity  with  phase  angle  0  <  f,  <  2rr,  and  Np  is  the  number  of  particles  in  the 
ensemble.  In  order  to  have  a  standard  comparison  between  various  sets  of  parameters  we  shall 
henceforth  consider  the  transverse  efficiency,  defined  by  %(;)  -  Ti(t)(y(ra)-\)/(utr0/c)2 . 
Thus,  r^ii)  is  the  efficiency  associated  with  the  transverse  beam  kinetic  energy  and  is  indepen¬ 
dent  of  v0j.  The  transverse  energy  efficiency  is  given  in  Fig.  (4)  as  a  function  of  time  for  the 
parameters  of  Fig.  (3).  The  net  energy  extracted  reaches  a  maximum  at  t  -  t,  -  165/<L;  the 
length  of  the  actual  cavity  would  be  L-v0lts.  Figure  (5)  shows  the  peak  transverse  efficiency  as 
a  function  of  «,  for  v^/c  -  0.3  and  various  values  of  the  frequency  mismatch  parameter 
Aw  —  w  -  ft  „/>(/■„).  From  these  curves  we  see  that  the  maximum  transverse  efficiency  is 
approximately  40%.  This  maximum  intrinsic  efficiency  of  approximately  40%  is  characteristic 
for  the  range  of  v0Jc  of  interest  in  the  maser.  We  note  that  in  our  treatment  the  steady  state 
value  of  «  is  taken  to  be  compatible  with  the  Q  of  the  cavity. 

III.  ENHANCEMENT  OF  MASER  EFFICIENCY 

A  qualitative  appreciation  of  the  nonlinear  particle  dynamics  can  be  gained  from  an  exam¬ 
ination  of  the  C-curves  and  phase  space  trajectories,  depicted  in  Figs.  2  and  3.  We  rote  that 
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since  the  particles  are  not  all  on  the  same  C-curves  they  reach  their  minimum  energies  at 
different  times.  Furthermore  the  minimum  energy  varies  widely  over  the  ensemble  of  parti¬ 
cles.  We  are  thus  led  to  the  conclusion  that  the  efficiency  can  be  dramatically  increased  by  first 
prebunching  the  particles  in  phase  space.  The  resulting  "macro-particle"  can  then  be  placed  on 
a  set  of  C-curves  for  which  virtually  100%  of  the  transverse  kinetic  energy  is  extracted. 


To  see  in  detail  how  this  can  be  accomplished  we  must  examine  the  basic  orbit  equations 
in  (7).  We  consider  the  linear  regime  of  these  equations  by  expanding  rand  £  in  powers  of  the 
normalized  electric  field  «.  For  e  sufficiently  small  we  can  write  r  =  r'0)  +  £r("  and 
f  =  £(01  +  e£,u  where  r0> ,  rU),  f‘0>.  are  independent  of  e.  The  solution  of  (7)  in  the 
linear  regime  is  r<0‘  -  rf),  £,m  =  -  A  ml  and 


—  (sin(f0  -  Aw/)  -  sin  f0) , 
2Aw 


f 


91  I 


2Awr0 


1  - 


fin  w2  ra2 


Aw  c2 


cu'  rf 

(costto  -  A w/)  -  cos  £,,)  +  n()/  — —  sin  £0 

r 


( 10a. b) 


where  |ertl)|  «  r{0) ,  and  |«£U)|  «  £‘01.  The  standard  deviations  of  rand  £  are  measures  of 
the  spreads  of  the  particle  ensemble  in  phase  space  and  are  defined  respectively  as 
o-,  —  (<r2>  -  <r>2)l/2  and  rri  —  (<£2>  -  <£>:>1,:  where  <(...)>  —  < 2-rr ‘ 

i*  2t r 

J </£„(...).  In  the  linear  regime  we  find 
<j,  “  a(l  —  cos  Ao >t)]  2 /■). 


rr2  . 

~3  2“ 


w2r02 


n0f  + 


1  - 


a0  w‘  rQ 


sin  Aw/i 


I  2 


<  1  la,b> 


where  a  -  «c/(2r0Aw)  will  be  denoted  as  the  "bunching  parameter."  The  mean  values  of  rand 


f  are  <r>  —  r0  and  <£>  —  rr  -  Aw/.  These  linearized  equations  are  valid  if  a  «  1  From  ■ 
Eqs.  (11)  we  note  that  a,  is  bounded  while  etc  decreases  secularly  in  time  Therefore,  if 
a  «  1,  the  particle  distribution  will  bunch  into  a  macro-particle  with  phase  space  dimensions 
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\r  ~  'J2ar0, 

=  (ir2/3  -  2a(<u/,0/c)2fl0i) l/2,  (12a, b) 

Note  that  the  expression  for  Af  is  invalid  for  times  such  that  fi0/  >  tr2(6a((or0/c):)_1.  The 
point  to  be  made  by  Eqs.  (12)  is  that  Ar  is  bounded  and  A£  initially  decreases  with  time  if  a  is 
sufficiently  small.  The  formation  and  dimensions  of  the  macro-particle  in  phase  space  is 
schematically  depicted  in  Fig.  6.  As  shown  in  Fig.  6,  all  the  particles  lie  on  open  trajectories 
and  the  formation  of  the  macro-particle  takes  place  over  many  cyclotron  periods,  with  only  a 
small  degree  of  spreading  in  r. 

We  can  tailor  the  bunching  parameter  such  as  to  achieve  minimum  <x(  when  the  macro¬ 
particle  enters  the  energy  loss  region  cos  f  >  0.  When  the  macro-particle  enters  this  region 
the  bunching  parameter  can  be  increased  such  as  to  deform  the  trajectories  of  the  macro¬ 
particle  onto  a  bundle  of  closed  C-curves  for  which  virtually  all  the  transverse  kinetic  energy 
can  be  extracted.  The  change  in  a  should  occur  over  an  axial  distance  roughly  equal  to  v0./Aw, 
in  order  to  ensure  that  the  macro  particle  remains  in  the  energy  loss  region.  The  bunching 
parameter  can  be  changed  as  a  function  of  axial  position  in  the  cavity  by  (i)  varying  K  as  a 
function  of  z  (for  example  by  axially  contouring  the  cavity  wall  radius)  and/or  (ii)  varying  the 
external  magnetic  field,  B0,  as  a  function  of  z.  From  the  definition  of  the  bunching  parameter, 
a  -  «c/(2r0Ao>),  we  see  that  these  methods  can  lead  to  a  large  fractional  change  in  a  with  only 
modest  changes  in  either  K  or  B0.  As  an  illustration  of  contouring  a  by  method  (ii)  we  choose 
the  magnetic  field  variation  depicted  in  Fig.  7.  The  initial  parameters  are  «  -  0.004,  -  0.3 

and  A(o/<L  -  0.08.  The  intial  value  of  the  bunching  parameter  is  a  -  0.077  (see  Fig.  7).  The 
particles  reach  their  maximum  degree  of  bunching  at  o> t  =  225,  just  before  entering  the  energy 
loss  region,  see  Fig.  8.  At  this  point  the  frequency  mismatch  is  decreased  gradually  to  the 
value  /a  -  0.02  at  &>/  •  260.  This  is  accomplished  by  increasing  the  external  magnetic  field 
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by  about  6%.  The  bunching  parameter  has  the  final  value  a  -  0.31  at  the  end  of  the  transition 
region  and  the  macro  particle  is  now  on  a  trapped  orbit.  The  maximum  transverse  kinetic 
energy  loss  occurs  at  d>r  -  400  resulting  in  a  transverse  efficiency  of  ■»  75%.  The  transverse 
efficiency  as  a  function  of  time  is  shown  in  Fig.  9.  As  shown  in  Fig.  5,  without  contouring  the 
applied  magnetic  field  the  transverse  efficiency  is  36%  for  these  parameters  with  constant 
A<u/£»  “  0.02. 

IV.  CONCLUSION 

We  have  shown  that  dramatic  increases  in  the  efficiency  of  the  maser  oscillator  can  be 
achieved  by  appropriately  tailoring  either  the  applied  axial  magnetic  or  the  cavity  field  (cavity 
wall  radius)  It  should  be  noted  that  this  approach  can  be  directly  applied,  with  appropriate 
modifications,  to  a  maser  amplifier  device.  We  have  illustrated  in  detail  an  example  for  which 
the  applied  magnetic  field  is  tailored.  The  requisite  field  variations  are  quite  modest  and  tech¬ 
nologically  straightforward.  No  attempt  has  been  made  to  obtain  an  optimum  field  profile,  as  it 
is  our  purpose  to  present  the  general  principle  of  efficiency  enhancement  and  show  that 
dramatic  increases  can  be  readily  achieved.  As  a  standard  of  comparison  we  have  considered 
only  the  transverse  efficiency,  which  is  independent  of  v0;.  In  principle,  the  energy  associ¬ 
ated  with  the  longitudinal  motion  may  be  largely  recovered  by  using  depressed  collectors,  so 
that  the  overall  system  efficiency  may  be  made  to  approach 
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Fig.  3  —  Phase  space  evoiulion  of  an  ensemble  of  16  particles  in  the  fields  of 
the  cavity  The  system  parameters  are  the  same  as  those  used  for  Fig.  1 2 )  a) 
initial  particle  distribution  bl  linear  regime  of  interaction,  c)  phase  bunching  re¬ 
gime.  d)  energy  extraction  regime  (Saturation) 
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Fig.  4  —  Transverse  energy  efficiency  as  a  function  of  lime 
(axial  position  within  cavity) 
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Fig.  5  —  Saturation  efficiences  as  a  function  of  the  normalized  cavity  field,  t.  for  various  values 
of  These  results  are  for  uniform  B0  and  K 
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Fig  6  —  Schematic  of  phase  space  evolution  illustrating  the  formation  of  a  macro-particie 
when  the  bunching  parameter,  a,  is  much  less  than  unity 


19 


STRANGLE  AND  SMITH 
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Fig.  7  —  Contoured  magnetic  field  profile  as  a  lunction  of  time  (axial  position  in  the  cavity)  in  the 
forward  wave  interaction  approximation.  The  final  parameters  correspond  to  the  uniform  field  sys¬ 
tem  of  the  previous  example  (Fig.  2-4) 


Fig.  8  —  Phase  space  evolution  for  (he  field  profile  shown  in  Kig.  7  a)  initial 
distribution  of  particles;  b)  formation  of  macro- pa ruc/e  ;usi  before  entering 
contoured  applied  magnetic  field  region;  c)  the  macro-panicle  after  ii  leaves  the 
contoured  applied  magnetic  field  region,  d)  particle  distribution  at  saturation 


Fig.  S  (Continued)  —  Phase  space  evolution  Tor  the  field  profile  shown  in  Fig. 
7  a)  initial  distribution  of  particles;  b)  formation  of  macro-particle  just  before 
entering  contoured  applied  magnetic  field  region;  c)  the  macro-particle  after  it 
leaves  (he  contoured  applied  magnetic  field  region,  d)  particle  distribution  at  sa¬ 
turation. 
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